Critical Current O-tt Transition in Designed Josephson Quantum Dot Junctions 
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We report on quantum dot based Josephson junctions designed specifically for measuring the 
supercurrent. From high- accuracy fitting of the current- voltage characteristics we determine the 
full magnitude of the supercurrent (critical current). Strong gate modulation of the critical current 
is observed through several consecutive Coulomb blockade oscillations. The critical current crosses 
zero close to, but not at, resonance due to the so-called O-tt transition in agreement with a simple 
theoretical model. 



Formation of a quantum dot (QD) between two su- 
perconductors, called a Josephson QD junction, enables 
study of supercurrent through a single energy level of an 
artificial atom (the QD). Supercurrent [Hill [3] through a 
QD m m El |7l E] is a dissipationless current which flows 
only when a constant phase difference across the QD 
is maintained [9l |10]. Uncontrolled fluctuations in the 
phase has so far prevented observation of the full mag- 
nitude of the supercurrent (the critical current). Here 
we analyze Josephson QD junctions, created in a sin- 
gle wall carbon nanotube, with a carefully designed on- 
chip circuit that controls phase fluctuations and thereby 
allows extraction of the critical current. The analysis 
reveals a reversal of the critical current for every elec- 
tron added to the QD, due to the so-called O-tt transition 
[IIl[Il[l3[ll[ia[Il[IIl[l8l[19], in agreement with a 
simple theoretical model. The realization of a well con- 
trolled nano-scale Josephson junction opens new routes 
toward superconducting quantum bits [20l |2TJ |22] . 

When a Josephson QD junction is placed in an elec- 
tromagnetic environment the phase difference across the 
QD (0) becomes a dynamical variable, which deter- 
mines the time- averaged current, /, and voltage drop, 
Vj = {{h / 2e) d(j) / dt) ^ across it. The /-Vj characteris- 
tics of the junction thus depends in a non-linear fash- 
ion on the phase dynamics induced by the electromag- 
netic environment. Inspired by Refs. [H [10] we utilize 
a designed external circuit in order to control the phase 
fluctuations which enables us to infer the true magni- 
tude of the critical current, /c, from the measurable crit- 
ical current /switching current, /m, by a theoretical fit- 
ting procedure. can significantly differ from 1^ as 
demonstrated previously for single wall carbon nanotube 
(SWCNT) based Josephson junctions [llSllZ]. 

To be able to design the external circuit we model our 
sample (Fig. [TJa) and (b)) by an extended resistive and 
capacitively shunted junction (RCSJ) model [4 , yielding 
the schematic circuit diagram shown in Fig. [ijc). The 
real Josephson QD junction is represented by an ideal 
Josephson junction with current-phase relation /j(0), in 
parallel with a junction capacitor Cj, and junction re- 
sistor i?j accounting for current carried by multiple An- 
dreev reflections. The resistor generally depends on 
Vj and gate voltage (V^ate) but for small enough voltage 
as used in our analysis Vj <C A/e, where A ~ 0.1 me V 
is the superconducting energy gap, we approximate i?j 



to depend only on Igate- Assuming a sinusoidal current- 
phase relation (/j((/>) = Ic sine/)) the dynamics of (j) in the 
circuit of Fig. [TJc) becomes equivalent to the damped 
motion of a fictitious particle in the so-called tilted wash- 
board potential [23 . The damping of the motion of this 
particle is characterized by the quality factor [H |10] 

^h[C{l^R/R,)^C,]/{2eQ 
^ RC^h/{2eI,Rj) ^ ^ 

given by the ratio of the local minimum oscillation fre- 
quency to the friction coefficient. Low Q implies higher 
probability of trapping the particle in a potential mini- 
mum resulting in a constant phase difference and, thus, 
the observation of the supercurrent. 

We have used Eq. ([T]) to design overdamped samples 
shown in Fig.JTJb). The superconducting electrodes (blue 
leads in Fig. [l[b)) are made with a small area to reduce 
the junction capacitance, Cj ~ 1 fF. Yellow leads, con- 
necting the superconducting electrodes to contact pads, 
are long thin normal metal wires designed to have a large 
resistance, measured to be ~ 1.5kQ. The contact pads 
are fabricated with a large area to increase their capaci- 
tance, C ~ 2pF. By inserting the above mentioned val- 
ues and Rj > h/e^ (see Fig. [sja)), into Eq. ([T]) we get a 
strong damping Q < 1 ioi 1^ > O.lnA which is also the 
smallest critical current we have been able to measure 
(see Fig.gb)). 

The measurements are performed in a ^He-^He dilu- 
tion fridge at a base electron temperature of 75 mK. In 
Fig. [ijd) we show a current biased four-probe measure- 
ment of an /-Vj curve at very low bias voltage V^d <C A/e. 
In this measurement, is the point where the curve has 
a large change in slope (as indicated in Fig. [TJd)), giving 
/m ~ 0.2 nA. For |/| < /m(< /c) the particle in the tilted 
washboard potential has a high probability after thermal 
activation out of one potential minimum to be subse- 
quently retrapped in the next potential minimum. The 
motion of the particle is therefore diffusive, leading to a 
small average phase velocity, i.e., a low but finite junction 
voltage. The branch in the /-Vj curve at |/| < /m, called 
the diffusive supercurrent branch, therefore has finite re- 
sistance. As the current is ramped further up / > /m 
the friction in the tilted washboard potential is no longer 
large enough to retrap the particle once it is activated 
out of a minimum. This leads to a high phase velocity 
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FIG. 1: Sample design and measurement scheme, (a) Schematic iUustration of the Superconductor(S)-Quantum dot(QD)- 
Superconductor(S) part of the device, made of a single wall carbon nanotube and a superconducting trilayer of 5nm Ti, 60 nm 
Al, and 5nm Ti, with transition temperature Tc = 0.75 K and energy gap A = 0.1 meV 5 . (b) Colored optical image of 
the device. Blue leads are the superconducting electrodes with a spacing of 300 nm, contacting a nanotube (represented by a 
dashed line). Yellow leads are long thin Cr/Au (15nm/15nm) leads connecting the superconducting leads to large- area contact 
pads, outside the image, (c) Schematic circuit diagram of the device with four-probe current-controlled measurement setup. 
Blue square: S-QD-S junction as depicted in (a). Yellow region: The resistors R are the yellow leads in (b), and the capacitors 
C are large-area contact pads, (d) Current (/) versus junction voltage (Vj) close to a resonance (see the corresponding arrow 
in Fig. |4jb)). Measurement (diamonds) are fitted with Eq. ^ (red line) yielding a critical current Ic = 1.3 nA and a junction 
resistance Rj = 90 kQ. 
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FIG. 2: Fitting procedure for extracting the critical 
current, (a) Theoretical current-phase relations calculated 
using a phenomenological "external exchange" model |17j , 
with ^ex = U/2 = 15A, r = Fl + Fr = llA, Fl = Fr, T = 
A/10, Iq = eA/h ^ 25 nA. The three curves are calculated 
very close to resonance at Ae = U/30, —U/60, —U/15 from 
top to bottom (see arrows in Fig. |4ja)), where As is the po- 
tential on the quantum dot measured from resonance. Even 
though this model is phenomenological it seems to capture 
correctly all qualitative features of the O-tt transition known 
from the numerical renormalization group 24 or quantum 
Monte Carlo [25 calculations. Diamonds in (b) are I-Vj 
curves calculated using the current-phase relations from (a) in 
the full theory [26], valid for a general current-phase relation 
and in the presence of thermal fluctuations. Eq. ([2| is fitted to 
the diamonds (solid lines) yielding very good fits in the whole 
gate voltage range apart from the closest vicinity of the O-tt 
transitions, which is below the experimental resolution. (See 
also the Supporting Information, Section S2). 



and hence the supercurrent is averaged out. 

In order to find Ic^ we fit the measured I-Vj curves 
to an overdamped (Q < 1) extended RCSJ model. For 



the overdamped Josephson junction we need to con- 
sider the classical dynamics only. In Ref. [27] the over- 
damped RCSJ model (without Rj) was calculated for si- 
nusoidal current-phase relation and was extended to gen- 
eral current-phase relations in Ref. [26^. High-resistance 
tunnel junctions have sinusoidal current-phase relations, 
but QD junctions with resistances comparable to the re- 
sistance quantum h/e'^ may have non-sinusoidal current- 
phase relation (see, e.g., Eq. (2) in Ref. [5]). Nevertheless, 
we approximate the current-phase relation by the simple 
sinusoidal form parameterized by /c, and we justify this 
approximation in Fig. [2] by comparing with a theoretical 
calculation where the full non-sinusoidal current-phase 
relation is included (see also the Supporting Information, 
Section S2). We generalize the theory of Ref. [27J to in- 
clude current carried via multiple Andreev reflections by 
the resistor Rj which is assumed to be much larger than 
the lead resistance R. The subsequent fitted values of Rj 
are consistent with this assumption, see Fig. ^c). Un- 
der these assumptions (/j(^) = Ic sin0 and Rj ^ R) the 
I-Vj curve is given parametrically via Vsd by the relations 



/(Fsd)=/clm 

^j(^sd) = ^sd- 



{I,h/2ekBT) 



{I,h/2ekBT) 



^j(^sd) 



(2) 



i^/(^sd), 



where /^(x) is the modified Bessel function of the com- 
plex order a, and rjiV^^) = hV^d/2eRkBT . Since all the 
parameters entering Eq. Q apart from 1^ and i?j are 
experimentally known we can use /-Vj curves to deter- 
mine both /c and i^j. The solid curve in Fig. [TJc) is a 
fit of Eq. ([2| to the measured data, yielding 1^ = 1.3 nA 
and Rj = 90 kQ. The theory seems to capture the ex- 
perimental measurement very well, which shows that the 
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FIG. 3: Supercurrent versus odd and even electron 
occupation, (a) Surface plot of the differential conduc- 
tance (dl/dVsd) versus bias voltage (14d) and gate voltage 
(V^gate). Coulomb blockade diamonds, indicated with black 
dotted lines, are alternating in size between large and small 
with a corresponding even and odd number of electrons lo- 
calized on the QD. (b) Surface plot of differential resistance 
(dVj/d/) versus applied current (/) and V^ate, in the same 
gate- volt age range as in (a), (c) Three I-Vj curves from (b) 
at indicated positions, the right and left graphs are shifted by 
20 /J.V for clarity. Circles are measured with odd occupation 
on the dot in the tt junction regime, squares at the O-tt tran- 
sition point, and triangles with even occupation on the dot in 
the junction regime. The solid lines are fits using Eq. S 
yielding critical currents of 0.65 nA(0.35 nA) for even(odd) 
electron occupation, and zero at the transition point. 



suppression of Ic = 1.3 n A into the diffusive supercur- 
rent branch with = 0.2 nA indeed is caused by ther- 
mal fluctuations. A total of five samples were fabricated 
with different contact resistances. They all showed I-Vj 
curves with qualitatively the same behavior as the sam- 
ple in Fig. [TJc) and showed excellent fits to Eq. ([2| (see 
Supporting Information, Section SI). 

In Fig. [3]^ a) we show a color-scale plot of the differential 
conductance {dl/dVsd) versus Vsd and l^ate- Coulomb 
blockade diamonds, indicated with black dotted lines, al- 
ternate in size between large and small as Igate is in- 
creased. This indicates that the QD has two-fold spin de- 
generacy of each discrete energy level, with an odd (even) 



number of electrons on the QD in the small (large) dia- 
mond. Coulomb repulsion energy U ^ 3meV and level 
spacing ~ 2.5 meV are extracted from this plot. The 
edges of the Coulomb diamonds are somewhat blurred 
due to the tunnel coupling between QD and leads. The 
tunnel coupling can vary from device to device and in 
the Supporting Information, Section SI, we show mea- 
surements on another device (fabricated the same way) 
exhibiting sharer edges of the diamonds due to a slightly 
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FIG. 4: Critical current O-tt transition compared to 
theory, (a) Solid line is |/c| versus potential on the quan- 
tum dot obtained as |/j(7r/2)| from [17] with the same pa- 
rameters as in Fig. [2] Crosses are the critical currents ob- 
tained by fitting Eq. ^ to the full theory as described in 
Fig. |2jb). Apart from very close to the transition point, the 
match between the two approaches is very good, (b) and 
(c) Experimental critical current |/c| and junction conduc- 
tance 1/Rj as extracted from fitting measured I-Vj curves to 
Eq. (U, versus the gate voltage V^ate- The error-bars in (b) 
at the plateau around Vgate = —2.56 V show the estimated 
precision of the fit. (d) Normal state differential conductance 
(B = 150 mT) at zero bias, in the same gate- voltage range 
as in (b) and (c). The dotted vertical lines indicate the po- 
sitions of charge-degeneracy resonances. Note the good cor- 
respondence between the phenomenological model in (a) and 
the measurement in (b), both the observed magnitude of the 
critical current as well as the shift toward odd occupation of 
all the O-TT transition points are in agreement with the theory. 
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lower tunnel coupling. At low bias voltage two paral- 
lel conductance ridges (indicated with black arrows) are 
seen, reflecting the peak in density of states of the super- 
conductors at Vsd = ±2A/e, yielding A ~ 0.1 meV. The 
energy broadening F = Fl + Fr, where ^l{r) is the cou- 
pling to the left (right) lead, and the asymmetry Fl/Fr 
of each Coulomb oscillation are determined in the normal 
state {B = 150 mT) by fltting the even valley part of the 
Coulomb oscillation peaks. Fig. 4(d), to a Lorentzian (see 
Supporting Information, Section S3). This gives approx- 
imately constant values of F ~ 1.1 me V and asymmetry 
parameters Tl/Tr in the range 1 to 4. These values are 
used in the theoretical plots in Fig. |2] and Fig. [4]^a). 

In Fig. [sj^b) we plot the differential resistance dVj/d/ 
as a function of / and Igate- Three representative I-Vj 
curves selected from Fig. IsFb) at indicated positions are 
shown in Fig. [sjc). The diffusive supercurrent branch 
is observed at most gate- voltages also at the resonances, 
but not in a narrow gate-region close to each resonance 
(open square in Fig.jsj^c) and Fig.|4jb)). Full gate- voltage 
dependence of /c, and Rj obtained by fltting the I-Vj 
curves from Fig. [sj^b) to Eq. ([2| is shown in Fig. ^h) 
and (c), together with the normal-state {B = 150 mT) 
zero-bias conductance in Fig.[4]^d). We observe that Ic os- 
cillates in accordance with the number of electrons (odd 
or even) localized on the QD. 1^ furthermore exhibits 
sharp dips to zero for every electron added to the QD, 
signifying a reversal of the sign of Ic due to a O-tt tran- 
sition of the current-phase relation (see Fig. [2| . All the 
transition points are seen to systematically shift toward 
odd occupation on the QD. In Fig. |4]^a) we have used 
the phenomenological model of Ref. [17 with the same 
parameters as used in Fig. |2] to calculate the full gate- 
dependence of the critical current (see Supporting Infor- 
mation, Section S2). Even though the QD in reality has 
an unpolarized magnetic moment in the odd valleys, the 
theory, which is based on a polarized magnetization, does 
capture two important consequences of the magnetic mo- 
ment: (i) the magnitude of the critical current and thus 
its strong suppression as compared to a non-interacting 
resonant level (/q = eA/h :^ 25 nA), which is a result of 
the transport of Cooper-pairs through a strongly corre- 
lated electronic system, (ii) the slight shift of the tran- 
sition points toward odd occupation. 

We end by noting that Josephson junctions fabricated 
in Al-AlOx-Al are already used in superconducting quan- 
tum bits (qubits) [20l [2TJ [22] . The realization of a well- 
controlled Josephson QD junction in a SWCNT, opens 
the route towards new superconducting qubits where bet- 
ter control of the Josephson junctions can be achieved. 
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To demonstrate reproducibility we show similar measurements on 
another device, (ii) Justification of the fitting procedure and ex- 
planation of the simple model prediction for the supercurrent. (iii) 
Estimating the asymmetry and level broadening of the device from 
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